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Abstract 

We perform an explicit two-loop calculation of the dilatation operator acting on single trace 
Wilson operators built from holomorphic scalar fields and an arbitrary number of covariant 
derivatives in M = 2 and M = 4 supersymmetric Yang-Mills theories. We demonstrate that 
its eigenspectrum exhibits double degeneracy of opposite parity eigenstates which suggests that 
the two-loop dilatation operator is integrable. Moreover, the two-loop anomalous dimensions 
in the two theories differ from each other by an overall normalization factor indicating that 
the phenomenon is not sensitive to the presence of the conformal symmetry. Relying on these 
findings, we try to uncover integrable structures behind the two-loop dilatation operator using 
the method of the Baxter Q— operator. We propose a deformed Baxter equation which exactly 
encodes the spectrum of two-loop anomalous dimensions and argue that it correctly incorporates 
a peculiar feature of conformal scalar operators - the conformal SL(2) spin of such operators is 
modified in higher loops by an amount proportional to their anomalous dimension. 
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1 Introduction 



The analysis of one-loop anomalous dimensions of composite high-twist operators in QCD re- 
vealed that its dilatation operator possesses nontrivial integrability symmetry. It was observed 
in Refs. [U El E] that the one-loop mixing matrix for the so-called maximal-helicity quasipartonic 
|4J Wilson operators can be mapped in the multi-color limit into the Hamiltonian of the SL(2; R) 
Heisenberg spin chain and its eigenspectrum can be computed exactly with the help of Bethe 
Ansatz. The length of the chain is determined by the number of elementary fields in the Wilson 
operator and the spin operators on its sites are defined by the generators of the collinear subgroup 
of the conformal group SO(4, 2). Let us emphasize that similar structures were discovered earlier 
in the Regge limit of QCD [5l Ej. Integrability imposes a very nontrivial analytic structure on the 
anomalous dimensions, which is reflected, in particular, in pairing of opposite-parity eigenstates 
in the spectrum. 

Integrability observed in QCD anomalous dimensions at one-loop order is a generic phe- 
nomenon of four- dimensional Yang-Mills theories and it is ultimately related to the presence of 
massless spin-one gauge bosons in the particle spectrum. It is also present, as was found in 
Refs. [TJ, [HI EJ , in maximally supersymmetric Yang-Mills (SYM) theory. Theories with less super- 
symmetries also inherit integrability although the number of integrable sectors strongly depends 
on the particle content of the models and is enhanced for theories with more supercharges |10j . 
In this regard the maximally supersymmetric M = 4 Yang-Mills theory currently occupies a 
distinguished niche in light of the gauge/string duality [TT] which allows one to establish the 
correspondence between the anomalous dimensions of composite operators in M = 4 theory and 
energies of string excitations on the AdSsxS 5 background p2J EH] • It was recently shown that 
classical string sigma models with anti-de Sitter space as a factor of the target space possess an 
infinite set of integrals of motion and therefore are integrable [T5| [16l [T7j. On gauge theory side, 
this suggests that the dilatation operator for Wilson operators carrying large quantum numbers 
should be integrable in the M = 4 theory in the strong coupling regime. 

The question arises whether the one-loop integrability of the dilatation operator and integra- 
bility of the classical string sigma model is a manifestation of the same universal phenomenon 
at weak and strong coupling, respectively, and if so then whether the "perturbative" dilatation 
operator exhibits integrability order- by-order in the coupling constant. Since the range of inter- 
action in the spin chain increases with order in 't Hooft coupling constant A = g^ M N c / (87i 2 ), - 
being merely nearest-neighbor at one loop, then stretching to three adjacent neighbors at two 
loops, etc. - the spectrum of anomalous dimensions should be determined by yet unidenti- 
fied long-range integrable spin chain. Recent extensive perturbative studies indeed support this 
conjecture [TH [HI EQl [2U [^ [23l [2H [251 [2S]. 

In this paper, we perform an explicit two-loop calculation of the dilatation operator acting 
on single-trace Wilson operators built from holomorphic scalar fields and an arbitrary number of 
covariant derivatives in M = 2 and M = 4 supersymmetric Yang-Mills theories. The motivation 
behind doing it is two-fold: (i) to establish the anticipated integrability as well as its dependence 
on the number of supercharges and (ii) to unravel the underlying integrable long-range interac- 
tion. In our previous publications [23] we have addressed the same questions in the sector of 
three-particle gaugino operators with arbitrary number of derivatives in super Yang-Mills theo- 
ries with M = 1,2,4 supercharges. We have computed the spectra of their anomalous dimensions 
and have found that the energy of the states with zero quasimomentum is double degenerate in- 
dicating the existence of higher conserved charges in addition to the quadratic conformal Casimir 



[27l l28l [T], El [18]. In the present paper we continue the analysis of noncompact sectors closed 
under renormalization and extend our consideration to single-trace Wilson operators built from 
L scalar fields X and arbitrary number of covariant derivatives, 

O n (Q) =tr{(iD + riX(0)(iD + ) n *X(Q)...(iD + r-X(Q)} , (1.1) 

where n = (m,n 2 , . . . ,til) denotes a set of L nonnegative integers. Here n^D^ = D + = d + — 
ig[A+, ■} is the covariant derivative projected onto the light cone with the help of the light-like 
vector n M . This projection automatically selects the maximal spin component. In order to avoid 
mixing with operators built from gauginos and gluons, we choose all X's to be the same and to 
possess the maximal charge with respect to the internal i?-symmetry group, that is, X = in 
M = 2 theory and X = <fi\ + i<f) 2 , conventionally called Z, in M = 4 theory. The scalar operators 
(11.11) carry the Lorentz spin N = m + • ■ ■ + and the canonical dimension N + L and they can 
mix under renormalization with operators having the same N and L. 

The paper is organized as follows. In Section [2] we employ a well-developed QCD technique 
for perturbative computation of the dilatation operator in the momentum representation. In this 
representation, the dilatation operator can be realized as an integral operator acting on light- 
cone momenta of scalar fields. The mixing matrix in the basis of local Wilson operators is simply 
obtained by forming Mellin moments of its integral kernel and the explicit expressions can be 
found in Appendix. In Section [3] we discuss the spectrum of anomalous dimensions for scalar 
operators (II. ip and describe integrable structures behind the two-loop dilatation operator using 
the method of the Baxter Q— operator. Section H] contains concluding remarks. 



2 Two-loop noncompact dilatation operator 

A concise representation of the entire tower of local Wilson operators is achieved by means of non- 
local operators with elementary fields located at positions Zi on the light-ray, X(zin^) = X(zi) 

O(z) =ti{X(z 1 )[z 1 ,z 2 }X(z 2 )[z 2 ,z 3 }...X(z L )[z L ,z 1 ]} , (2.1) 

where z = (zi, z 2 , . . . , z L ) and the Wilson line [zj, Zj + i] = ig f* 3 ^ dzA + (z) is stretched between 
two fields to make the composite operators gauge invariant. The local Wilson operators (11.11) are 
deduced from O(z) by means of the Taylor expansion 

0W = f; WW ...H^OU(0). (2.2) 
' ml n 2 \ n L l 

m,n2,---,n L >0 

For our purposes yet it will be extremely useful to use a representation of the same operators in 
the reciprocal momentum space. It is given by the Fourier transform of O(z) with respect to the 
light-cone coordinates 

/oo j roc j 

-oo 27r J-oo 27T 

with u = (u±, u 2 , . . . , ul) being the vector of the light-cone momenta. These operators obey the 
renormalization group (Callan-Symanzik) equation 

(v^ + PM-jjjj) 5 H = / idv} L nu\v)0(v) , (2.4) 
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with the integration measure [(1v]l = dv\ . . . dviS(^2 k v k — J2 m u "i) and the dilatation operator 
in the momentum-space representation, Y(u\v), admitting perturbative expansion in 't Hooft 
coupling A = gl M N c /(Six 2 ) 

Y(u\v) = XY m {u\v) + X 2 Y (1 \u\v) + C(A 3 ) . (2.5) 

A detailed account on the technique used for the perturbative calculation of the kernels Y^(u\v) 
and Y^(u\v) can be found in our previous publication [23] . 

In this paper, we calculate the two-loop kernel Y(u\v) for nonlocal light-cone scalar operators 
( 12. ip in supersymmetric Yang-Mills theories with M = 2 and M = 4 supercharges. We shall 
employ the light-cone formalism of Ref. [10J which heavily relies on the light-cone gauge A+(x) = 
and which allows one to treat all supersymmetric theories in a unified fashion. Under this 
gauge condition, the Wilson lines in (12.11) reduce to unity [zfc,;Sfc+i] = 1 and the number of 
relevant Feynman diagrams decreases significantly. In the light-cone formalism, the Lagrangian 
of super Yang-Mills theory depends on physical components of elementary fields only since all 
non-propagating degrees of freedom can be integrated out. This allows one to switch from 
covariant spinor and vector fields to single-component fermionic and bosonic fields carrying a 
definite helicity and, then, introduce dimensional regularization inside Feynman integrals by 
continuing them to D = 4 — 2e dimensions without breaking supersymmetry of the underlying 
gauge theory. The resulting regularization procedure is equivalent to the dimensional reduction 
scheme. 

2.1 Evolution kernel 

Let us start with the one-loop kernel Y^(u\v). In the multi-color limit, it has a simple nearest 
neighbor structure 

vcvw=E|<Ui- r(o) ^-^)} n ( 2 - 6 ) 

k=l ^ ' i=i. 

with the periodicity condition L + 1 = 1 and the total momentum conservation J^ fc = Yl m Um 
absorbed into the integration measure in (12 .4p . Here, r^ ^ = |( jV — 4) is an additive constant 
and the two-particle kernel is defined as 

v i°fc+i = \®{ u k,v k )f s {u k ,v k )} { + k) + [G(u k+ i,v k+ i)f s (uk+i,v k+1 )]^ k+l) , (2.7) 

where the generalized step-function Q(u,v) specifies possible values of the momentum fractions 
and a notation was introduced for the decoration factor f s , 

Q(u,v) = 9(u)9(u-v)-9(-u)9(v-u), f s ( u ,v) = . (2.8) 

v — it 

The plus-distribution in (12.71) regularizes the end-point singularities of f s (u, v ) as u — v — » and 
is conventionally defined as 

= ^ - 5{u - v) [ dw^L , (2.9) 
v — u v — u / v — w 
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Figure 1: Different topologies of two- loop diagrams in the light-cone superspace formalism of Ref. 
[TO] involving three nearest superfields. The three-particle kernel (12.111) is obtained by projecting 
the superfields onto the scalar field component. Mirror symmetric diagrams are implied. The 
graphs with two-to-one and one-to-two particle transitions all vanish. 



for any test function r. 

We now turn to the two-loop contribution to (12. 5p . In the multi-color limit, the interaction 
can involve up to three nearest neighbor fields and the generic structure of the kernel Y^(u\v) 
reads 



¥ {1) {u\v) 



L 

E 

k=l 



V 



(i) 

k,k+l,k+2 



+ 5{u k -v k ) 



V 



(1) 

k+l,k+2 



n 



fe+l,fe+2 



(2.10) 

The three- and two-particle kernels 'Vk,k+i are found from the diagrams displayed in 

Figs. [T] and [2J respectively. Finally, receives contributions from two-loop renormalization of 
the external scalar legs. We will not dwell on the computation procedure since it was spelled out 
in details in Ref. |23j . 

The calculation of the Feynman graphs in Fig. [1] yields the three-particle kernel in M = 2 
and Af = 4 SYM theories as 



V 



(i) 

123 



©m/sK; v i)f s (u 3 , v 3 ) In 



(u 3 ) 



+ 



® ( 12 3 fs(Ul,Vi)f s (u 3j V 3 ) In 



vi(u 3 -v a - v 2 )u 2 _ 

{Vl + ^3 - ^l)(f2 +V3- U 3 ) 
U 2 V 2 



(uiu 3 ) 



(2.11) 



U\ <-> u 3 
Vl v 3 

In general, the other kernels Vj^L. x k+2 can be obtained from this expression through cyclic 
permutation of indices under the periodic boundary condition L + k = k. Here the generalized 
three-particle step-functions are defined following the conventions of Ref. [23] as 



e 
e 



(i) 

123 

(4) 
123 



Q(u h Vx) [Q(u 2 , Vi + V 2 - Mi) - 6(^2 + v 3 - u 3 , v 3 )} 
9(m 2 , v 2 + v 3 - u 3 )Q(v 2 + v 3 - u 3 , v 2 ) . 



(2.12) 
(2.13) 



For v 1,2,3 > and v\ + v 2 +v 3 = 1 they define the regions in (ui, u 3 )— plane (with u 2 = l — Ui — u 3 ) 
which produce non- vanishing contributions to 







(i) 

123 



< u± < vi 

V2 + V 3 < U 3 < 1 



Ui 







(1) 

123 



Vl < U\ < 1 — u 3 
v 3 <u 3 <l 



(2.14) 



4 




Figure 2: Topologies of two-loop Feynman graphs providing two-particle contribution to the 
dilatation operator. The blobs stand for one-loop correction to the corresponding propagators 
and three- and four-particle interaction vertices. Mirror symmetric diagrams have to be added. 



Remarkably, the scalar kernel (12. lip can be obtained from the three-particle gaugino kernel (see 
Eqs. (5.7)-(5.14) in |23j) by substituting the gaugino decorating factor f q (u,v) = u/(v(v — it)) 
with the scalar one f s (u,v). 

In Eq. (12.111) . in the first term the superscript (w 3 ) indicates that the plus-distribution is 
taken only with respect to the variable u% as shown above in Eq. (12. 9p . The second term in 
(12. lip involves the double plus-prescription defined as 



t(u, u') 



(v — u){v' — u') 



(uu 1 ) 



t(u, u') 



(v — u)(v' — vl) 
— 8{u — v) dw- 



+ 8(u — v)8(u — v) Idwdw' 



t(u, w 



(v — u)(v' — w') 



t(w, v!) 



w)(v' 



u' 



8{u' — v') I dw 1 



t(w, w') 



(v — w){v' — w') 



with t(u,u') being a test function. 

The two-particle kernel ~V^ k+1 receives contributions from the diagrams presented in Fig. |2] 
as well as from subtraction terms coming from the single and double plus-distributions. The 

12 



expression for looks like 



V 



(i) 



12 



M h - In 2 In — In 

6 4 v\ vi 



1-^ 



(«0 



+ 



Ul <-> U2 
V l J J J + L V l V 2 

(2-15) 

As before, all other can be obtained from by cyclically permuting the indices. Finally, 

one has to add a single-particle contribution due to two-loop renormalization of the scalar fields. 
The former can be read off from the renormalization constant of the gauginos in the dimensional 
reduction scheme discussed at length in Ref. [23J 



rW = (AT- 4) (2 -AO 



(2.16) 



so that = in J\f = 2 and M = 4 theories. This result has a simple interpretation (see 
Eq. fl2~2T]) below). 

It is interesting to notice that the three-particle kernel V^3 does not depend on the number 
of supercharges M. At the same time, examining the M— dependence of the two-loop kernel 
one finds the following remarkable relation between the two-loop dilatation operators in M = 2 
and Af = 4 theories, 

^M=2(u\v) = (1 + 2A) Y M=4 (u\v) . (2.17) 
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The same relation was previously found for the gaugino operators in Ref. [23]. It suggests that 
the dilatation operators in the two theories share the same properties to two loops at least, and 
their eigenvalues are, obviously, related via Eq. (I2.17p . 

Being combined together, Eqs. (12.51) . (I2.10p . (12. lip . (12.151) and (12.161) provide explicit expres- 
sions for the two-loop dilatation operator for the scalar operators (12. ip in M = 2 and M = 4 
SYM theories in the multi-color limit. In section [HI we will diagonalize the dilatation operator 
(12.51) and determine the spectrum of anomalous dimensions of scalar operators. To accomplish 
this goal however, we will first reconstruct the mixing matrix for the operators (11.11) . 



as 



2.2 Mixing matrix 

It follows from (I2.2p and (12.31) that the local Wilson operators (11.11) are related to the moments 
of the operators in the momentum representation 0(v) 

0,(0) = J dvtv™ 1 ...J dv L v n L L 0(v) . (2.18) 

Together with (12.41) this implies that the operators O n (0) obey the renormalization group equa- 
tion 

with m = (mi, . . . , mi) such that m 8 > and J2 k = J2 k m^. The mixing matrix is introduced 
I [du] L uTuf ...u n L L Y(u\v) = -J2 A™(A)<^ 2 m2 ■ ■ . . (2.20) 

It defines the representation of the evolution kernel Y(u\v) in the space spanned by the poly- 
nomials v i 1 v 2 2 ■ ■ ■ v ™ L . Substituting the evolution kernel in (I2.20p by its two- loop expression 
(12.51) . we obtain the corresponding mixing matrix in the multi-color limit. In particular, for 
n\ — . . . — til — one finds that the terms in (12. 5p containing '+' and '++' distributions provide 
vanishing contribution to the left-hand side of (12.201) leading to 

A°(A) = L [Ar<°> + A 2 r^) + 0(A 3 )] = | [X(Af - 4) + A 2 (A/" - 4) (2 - AT) + 0{\ 3 )} (2.21) 

According to (I2.19P and (12.201) . Aq(A) defines the anomalous dimension of the local operator 
O(0) = tr [X L (0)] in M = 2 and M = 4 theories. These operators are known to be protected [29] 
and, therefore, their anomalous dimension is equal, in our notations, to L/3^(g)/2. This implies 
that (12.211) is exact to all loops. 

The perturbative expansion of the mixing matrix A^(A) is similar to that of the evolution 
kernel (12.51) and reads in the multi-color limit 

A^(A) = AA (0) ^ + A 2 A (1) ^ + 0(X 2 ) (2.22) 
with A = g 2 N c /(M 2 ). The matrices A^ and A^ are given by the moments (12.201) of the kernels 
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V(°)(v|i>), Eq. (EE}, and VW(w|v), Eq. (Efl}, respectively, 



a (l)m __ |,(l)m t m H1 m H2 . (l)m fc m fc+1 ™ t+2 p(l)rm fc rm fc+1 rm i+2 1 TT rrry 

^\ n k n k+1 n k+2 " rJl "fcWfc+i "n fc+ 2 ~i~ - 1 - "n fc "n fc+1 "n fe+2 J XI ™J ' 
fc=l i=t,iitk, 

k+X,k+2 

Here 5™ is the Kronecker symbol and periodic boundary conditions L + k = k are implied. We 
also introduced a notation for the moments of the two- and three-particle kernels. The explicit 
expressions for the corresponding matrices A I™ 1 " 12 , A^" 11 " 12 and A^™ 1 ™ 2 " 13 are rather lengthy 
and can be found in Appendix. 



3 Eigenspectrum and integrability 

To solve the evolution equation ( 12. 19ft . one examines the eigenproblem for the mixing matrix 

(EL22D 

£A£(A)*(n)= 7 (A)*(m), (3.1) 

where the sum runs over n = (ni, . . . , til) such that rij > and J2k n k = J2k m k- We remind 
that the (3— function vanishes in the Af = 4 SYM theory /3^=4(^) = 0, while its exact value 
in the Af = 2 SYM theory is given by the one- loop expression /3jv =2 (g) = —2 A. Then, it 
follows from Eq. (12. 19j) that in Af = 4 theory, the conformal operators C con f (0) = J2 n ^( n )^n(0) 
have an autonomous scale dependence with 7(A) being the corresponding anomalous dimensions. 
In Af = 2 theory the scale dependence of the operators C con f(0) is more involved due to the 
dependence of the eigenstates \l/(n) on the running coupling constant but one still refers to 7(A) 
as an anomalous dimension. 

Having the explicit expression for the two-loop mixing matrix (12.221) at our disposal, we can 
solve the spectral problem (13. ip for arbitrary lengths L and look for manifestation of symmetries 
of the dilatation operator in its spectrum. The evolution kernel V(tt|v) preserves the total 
momentum P = ^2 n u n = ^2k v k an d its eigenvalues 7(A) do not depend on P. This allows 
one to simplify the analysis by going over to the forward limit, P = 0, see Ref. [23]. Still, 
the anomalous dimensions 7(A) depend on the total number of derivatives N = n i which 
is one of the integrals of motion for the Schrodinger-like equation (13.11) . Another conserved 
charge follows from the invariance of the mixing matrix (|2.22[) under discrete cyclic P and mirror 
transformations M defined as 

P*(ni,n 2 ,...,n L ) = *(n 2 , n 3 , . . . , n x ) , (3.2) 
M*(m,n 2 , . . . ,n L ) = ty(n L ,n L -i, ■ ■ ■ ,ni) • 

Since these two operators do not commute with each other, the solutions to (13.11) can be classified 
according to the eigenvalues of only one of them, say P 

P*(m, n 2 , ...,n L ) = e ie ^(n l7 n 2 , . . . , n L ) , (3.3) 
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(b) 



10 15 20 



(c) 



Figure 3: The spectrum of anomalous dimensions of L = 3 operators (see text) at one loop (a), 
and two loops for Af = 2 (b) and M = 4 (c) super Yang-Mills theories. 



where the quasimomentum 9 takes L distinct values, 9 = with n = 0,1, . . . , L — 1. Making 
use of the relation PMP = M, one immediately finds from (13.11) that for 9^0 the eigenstates 
\P(n) and M\E'(n) have the same "energy" 7(A) and opposite values of the quasimomentum, 9 
and —9, respectively. Thus, the solutions to (13. ip with nonzero quasimomentum are necessarily 
double degenerate. For the eigenstates with 9 = the discrete symmetry alone does not imply 
any degeneracy. We recall that the eigenstates \l/(n) determine the form of conformal operators 
Cconf(O) = J2 n ^(n)On(0) with the basis operators O n {0) defined in (II. ip . Since the latter 
operators are cyclically symmetric with respect to n, the eigenstates \I/(n) should possess the 
same symmetry, that is, they ought to have zero quasimomentum 9 = 0. This leads to an 
additional selection rule for solutions to Eq. (13. ip . We would like to stress that the mixing 
matrix (I2.22p possesses eigenvalues with both quasimomenta 9 7^ and 9 = but only the latter 
define eigenvalues of the dilatation operator in gauge theory. 

To two-loop accuracy, the anomalous dimensions in the multi-color limit have the perturbative 
expansion 

7 (A) = A e (0) + A 2 + C(A 3 ) , (3.4) 

with e^ ) and e^- 1 ' being functions of the length of the operator L, the total number of derivatives 
N and some other quantum numbers, yet to be determined. To find the explicit form of these 
functions one has to diagonalize the mixing matrix for various L and N. In particular, using the 
expression for the two- loop mixing matrix, Eqs. (I2.22p and (12.241) . and solving the eigenproblem 
(13.11) for scalar operators of length L = 3 and the total number of derivatives < iV < 20, we 
calculated the values of and with 9 = in the Af = 2 and Af = 4 SYM theories. We 
summarized our results in Fig. [31 We found that in both theories all eigenvalues (13 .4p (except of 
a single lowest eigenvalue for each even N) are double degenerate to two loops . We would like 
to stress that the two-loop evolution kernel in the Af = 2 theory contains conformal symmetry 
breaking terms proportional to the (3— function. In the same fashion as in Ref. [23], the fact that 
the degeneracy is present in the Af = 2 theory suggests that the phenomenon is not directly tied 
to the conformal symmetry. 
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3.1 One loop 

We recall that for eigenstates with vanishing quasimomentum 6 = 0, the discrete symmetry 
(13.21) is not sufficient to warrant the double degeneracy of the spectrum. To one-loop order, this 
property immediately follows from hidden integrability of the one-loop dilatation operator [271 [U 
El [18] . Namely, the one-loop mixing matrix 

A (o) 

can be mapped into a Hamiltonian of the XXX 
Heisenberg spin chain of length L and spin operators in all sites being generators of (infinite- 
dimensional) discrete series representation of the 5L(2;R) group of spin j = |. On the gauge 
theory side, this group emerges as the so-called collinear subgroup of the full (super) conformal 
group. Its representation is determined by the conformal spin of scalar fields entering the Wilson 
operators (II. ip . It is important that for the one-loop dilatation operator, the conformal spin of 
scalar field can be substituted by its value at order A , that is, by its classical value j = ~. As 
we will argue below the situation becomes more complex starting from two loops. 

Thanks to integrability, the spectral problem (13.11) can be solved exactly to one-loop order 
using the Quantum Inverse Scattering Method [30J. This method allows us to identify the 
complete set of conserved charges qf \ ■ ■ ■ , q® whose eigenvalues parameterize solutions to (13. ip . 
For scalar operators of length L, they can be obtained as coefficients in the expansion of the 
auxiliary transfer matrix in powers of the spectral parameter u 



t^\u) = tr {L l («)Lh(m) . . .Li(u)} = u L |2 + q^vT* + q [ i"u~' d + ■■■ + qf'u^l . (3.5) 

Here Lfc(u) = t ■ u + ia ■ Sk is the standard Lax operator for the XXX Heisenberg spin chain 
and a a are Pauli matrices. It is given by a 2 x 2 matrix whose entries are spin operators S% in 
k th site. In gauge theory, the latter are generators of the collinear SL(2) subgroup acting on k th 
scalar field in (II. ip . The "lowest" integral of motion q^ } is related to the total conformal spin of 
the Wilson operator, J = N + L/2, 

qf = -{N+\L){N+\L-l)-\L. (3.6) 

The explicit form of higher conserved charges gj^ 2 can be found in [27]. In particular, they 
possess a definite parity with respect to the discrete transformations (13.21) . [P, q^] = {M, g^fc+i) = 
[M, g|fc ] — 0. This implies that the one-loop "energy" e(°) as a function of the conserved charges 
q% satisfies the relation 



(°)„,-2 i J°)„,-3 , , „(o)„,-£ 



e 



(°> (qf\ qf\..., q f ) = e <°> ( - qf\ qf\ . . . , (-l)^f ) , (3.7) 



and, therefore, all eigenstates including those with zero quasimomentum states are double de- 
generate provided that q^ +1 ^ 0. By explicit diagonalization of the dilatation operator in M = 2 
and Af = 4 SYM, evaluated in the previous section, we found that its eigenspectrum (13.41) is 
double degenerate not only at one loop but also in two-loop order. Together with (13 .7p this 
suggests that the all-loop anomalous dimension (13 .4p is a function of the conserved charges gfc(A) 

q k (X) = qf + Xq ( k 1) + 0(X") } (3.8) 

such that it verifies the same parity relation (13.71) . We shall argue in Sect. 13.21 that this is indeed 
the case to two loops, at least. 
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The eigenspectrum of the one-loop anomalous dimension can be easily found using the 
method of the Baxter Q— operator for the SL(2; R) Heisenberg spin chains. The method relies 
on the existence of a commuting family of operators Q(w) which depend on an arbitrary spectral 
parameter u and, in addition, commute with the Hamiltonian (i.e., the one-loop dilatation opera- 
tor). The one- loop anomalous dimensions and the corresponding quasimomenta 9, Eq. (13. 3p . 
can be expressed in terms of the eigenvalues of Baxter operator that we shall denote as Q^°\u) 



e 



i (W 0, (i)) - i M (0) H )) , ^ = oWpiy ' ( 3 - 9 ) 



For the SL(2;M) magnet of spin j and length L the function Q^(u) satisfies the second order 
finite-difference equation [31], the so-called HQ" or Baxter equation 

Af (u)Q {0) (u + + A {0 \u)Q {0) (u -i) = tf\u) Q {0) (u) , (3.10) 

where the "dressing factors" A±\u) = (u±ij) L depend on the spin j and t£ (u) is a polynomial 
in u of degree L defined in (I3.5p . We recall that for the scalar operators (II. ip . to one-loop order, 
the spin j is given by the conformal spin of scalar field at A order, i.e., 3 = \- The Baxter 
equation 1 13. 10p alone does not specify Q^(u) uniquely and it has to be supplemented by an 
additional condition that Q(°'(u) should be a polynomial in u of degree N > [27J. Then, 
Q(°'(u) can be parameterized (modulo an overall normalization) by its roots 

N 

Q {0 \u) = H(u-uf ) ). (3.11) 

k=l 

It is known that for the SX(2;R) spin chain the roots take real values only [27]. In gauge 
theory, the nonnegative integer iV coincides with the total number of derivatives in (11.11) . Being 
combined together, Eqs. (13.101) and (13.111) uniquely define Q^°\u) and allow one to calculate the 
one-loop anomalous dimensions (13.91) (see Fig. [3^,) and determine the corresponding values of the 
conserved charges qf 1 entering (13.51) . 

The method of the Baxter Q— operator is equivalent to the Bethe Ansatz. Indeed, substituting 
u = Ufc into both sides of (13 . 1 Of) . one finds that u£ satisfy the Bethe-root equations. In the 
same manner, substitution of (13.111) into (13.91) leads to the well-known expressions for energy it 
terms of Bethe roots. In the next section, we will discuss a "deformation" of the one-loop Baxter 
equation (13.101) in order to accommodate the two-loop corrections to the anomalous dimensions 
that were computed in Sect. [2j The reason why we prefer to deal with the Baxter Q— operator is 
that its eigenvalues have a direct physical meaning which should be preserved to all loops - for 
real u the function Q(u) determines the wave function of the spin chain in separated ("collective") 
variables [321 [33] (see also discussion in Sect. 4). Defined in this way, Q{u) should oscillate on 
the real u— axis and the number of its nodes should be equal to the excitation number N. 



3.2 Two loops and beyond 

The double degeneracy of the spectrum of two- loop anomalous dimension established in Sect. [2] 
combined with the exact integrability of the one-loop spectrum suggests to generalize the Baxter 
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equation to higher loops. Going over to two loops, we expect that the Bethe roots will be 
corrected by perturbative corrections 

Mfc (A) = 4 0) + A4 1) + O(A 2 ) (3.12) 

and Uk will verify "modified" Bethe equations. For the scalar operators under consideration, such 
equations have been conjectured in [211 E3! 



X +\ L JL x - - x + 1 - A/ 2x+xT 

-M =114 w 9 - + . ( 3 - 13 ) 



where x^ = x(uk ± f ) and the deformed spectral parameter x = x(u) is defined as [35J 



x{u) = \u 1 + y/1 - 2\/u 2 . (3.14) 

Then, the anomalous dimension and the corresponding quasimomentum are determined in terms 
of x^ parameters as follows [211 EI] 

N / • • \ N + 

This relation coincides with (13 . 91) to one loop and it is believed that it should reproduce the 
anomalous dimensions of scalar operators (11. ip in the M = 4 SYM theory to three loops at least. 

As was already mentioned above, the one- loop Bethe equations for the parameters u% are 
equivalent to the Baxter equation (I3.10P for the polynomial Q(u), Eq. (13.111) . We assume that 
the same relation between the Bethe Ansatz and the Baxter equation also holds in higher loops 
and introduce into consideration a polynomial with roots given by parameters U)~, Eq. (I3.12p 

N 

Q{u) = l[(u - u k (X)) = Q(°\u) + A QW(u) + \ 2 Q {2 \u) + .... (3.16) 

k=l 

Here the leading term Q (°>(u) coincides with the solution to the Baxter equation (I3.10P and is 
given by a polynomial of degree N with real roots uf\ By construction, the subleading terms 
Q( n \u) do not depend on the coupling constant and are given by polynomials of degree N — 1, 
that is, Q {n) (u) ~ u L ~ l for n > 1. 

It turns out that the function Q[u) defined in (13.161) obeys a second-order finite difference 
equation very similar to the Baxter equation (13.101) 

A+(x(u+i))Q(u + i)+A_(x(u-i))Q(u-i)=t L (x(u)) Q(u) , (3.17) 



where the dressing factors A a (x) (with a — ±) satisfy the condition A_(x) = A + (x*) and read 
to three-loop accuracy 

A a (x)=x L A^(x), (3.18) 
with the function x = x(u) defined in (I3.14p and 



A^Hx) = exp f~ (lng(fa))' - ^ (inQ^a))" + x (inQ(fa)) 



+ G(A 3 ) . (3.19) 
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To lowest order in A one has x(u) = u + 0(A) with A ( ± n \x) = 1 + 0(A) leading to A±(u) = 
(u±|) L . It is straightforward to verify that the roots of Q(u) verify the modified Bethe equations 
( I3.13P to three- loop accuracy [21]. The apparently unusual feature of (13.171) and (13 . 18[) compared 
to (13. 10p is that the dressing factors A + (u) and A_(w) depend on derivatives of the function 
Q(u) evaluated at u — | and u — — |, respectively. We shall elucidate the origin of this property 
in a moment. 

The auxiliary transfer matrix £l(w) entering the right-hand side of (13. IT)) is a generating 
function for the conserved charges (13.8)) 



t L (x) 



X 



(3.20) 



n>2 



Its perturbative expansion starts with (I3.5P and includes higher-loop perturbative corrections to 
the integrals of motion. In distinction with (|3.5j) . the series in the right-hand side of (13. 20)) does 
not truncate, thus, reflecting the asymptotic character of the Baxter equation ( 13. 1 Tf) . Notice the 
series in the right-hand side of ( 13.201) does not involve ~ x~ x term. To see this, it suffices to sub- 
stitute (13.161) and (13.20)) into (13. IT)) . expand its both sides at large u and match the coefficients 
in front of powers of u. In this manner, one deduces from Eqs. (13 . 1 T)) that the 'lowest' charge 
92(A) equals 

<fe(A) = -{N + \L + | 7 (A)) (N + \L + | T (A) - l) - \L , (3.21) 



where 7(A) is given by the sum of two functions 7 CT (A) (with a 
asymptotic behaviour of (x) at large x 



±) parameterizing leading 



7 (A) = 7+ (A) - 7 -(A) , A^\x) = exp (z 7t7 (A) x~ l + 0(oT 2 )) . 



(3.22) 



They satisfy the relations 7+ (A) = —(7- (A))* and 7<r(A) = 0(A) which ensure that 7(A) takes 
real values and vanishes for A — > 0. As a result, the charge 92(A), Eq. (13.211) . also takes real 
values and approaches its lowest order value (1 3 . 6 p for A = 0. The reason we used in (I3.22p the 
same notation as in (I3.15P is that 7(A) defines the multi-loop anomalous dimension. Combining 
together (13TT5D and (ET22]I we get 



7(A) 



X— —— — — 

du 4 du 3 48 du 5 



(ilnQ(w)) 



u=i/2 



A£ (o) + \ 2 £ ( 1 ) + \z £ m + . . . . (3.23) 



-i/2 



To leading order in A, it coincides with (13.91) and matches (I3.15P up to three loops. One can apply 
Eqs. ( 13.171) and (I3.23P to determine the spectrum of anomalous dimensions of scalar operators in 
M = 4 SYM theory, Eq. (11.11) . of arbitrary length L > 2 and total number of derivatives N > 0. 

Solving the Baxter equation (|3.17j) . we can determine the function Q(u) up to three- loop 
order, Eq. (I3.16p . as well as the conserved charges 9fc(A), Eq. (13.81) . In particular, gfc(A)'s take 
real values and Q(u) is a real function of u. The fact that the leading function Q(°\u) has 
only real roots ensures that the three-loop Bethe roots Mfc(A) are also real. The corresponding 
quasimomentum is given by the leading order relation (13. 9p involving the Q^°>— function and it 
is protected from perturbative correction in A. 

To elucidate the physical meaning of the factor A^ ron) (x), Eqs. (13181) and (13191) . it is in- 
structive to compare the SL(2) Baxter equation (13.171) with a similar equation describing the 
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multi-loop anomalous dimension of Wilson operators in the 577(2) sector in M = 4 SYM the- 
ory [35], i.e., single-trace operators of canonical dimension L built from iV holomorphic scalars 
X = 03 + an d L — N fields Z = (f>i + z</> 2 . In that case, the SU(2) Baxter equation takes 
the same form as (13 . 1 71) with the only difference that the "dressing" factors are merely given 
by A± {2) (u) = (x(u =F f )) L and do not involve additional factors similar to Aa en \x). A natural 
question arises what is the reason for such difference? We recall that the anomalous dimensions 
( 13.41) are eigenvalues of the dilatation operator which, in its turn, is one of the generators of the 
(super) conformal group of the underlying gauge theory. In the SU(2) sector, the mixing occurs 
between Wilson operators carrying the same canonical dimension L and the isotopic charge N. 
Then, in the multi-color limit, the dilatation operator in this sector can be mapped into the spin 
chain with the spin operators being the generators of the 577(2) subgroup of the full i?-symmetry 
group. This should be compared with the SL{2) sector in which case the dilatation operator is 
mapped into the spin chain in such a way that the spin operators are generators of the collinear 
57(2) subgroup and the dilatation operator is one of these generators! 

To lowest order in A the recursion works as follows. The one-loop dilatation operator is 
identified as a Hamiltonian of the SL(2) spin chain with the spin operators being the generators 
of the collinear subgroup in gauge theory to zero-loop order0 To this order, the generators of the 
collinear subgroup depend on the classical value of the conformal spin of the scalar field j = ~ 
and, as a consequence, the one-loop dilatation operator coincides with the 57(2) spin chain of 
spin j = Going over to higher orders one expects that the dilatation operator to n th loop is 
given by the SL(2) spin chain with spins corrected by perturbative corrections to the dilatation 
operator to (n — l) st — loop accuracy. This property finds its manifestation in the structural form 
of the Baxter equation ( 13.171) . 

Indeed, the Baxter equation (13.171) involves the dressing factors Aa en \x), Eq. (13.191) . which 
depend on the Q— function that satisfies the Baxter equation itself. Replacing Q{u) in (13.171) by 
its perturbative expansion ( 13.161) and expanding A& en \x) in powers of A it is easy to see that 
Aa Bn \x) induces corrections to the Baxter equation for Q^ n \u) involving lowest-order functions 
QW(m) with < k < n. In particular, for n = 1 the one-loop corrections to the conserved charge 
#2 (A) depend on the Q^>— function, Eq. ( I3.2ip . We remind that the charge qf\ Eq. ( 13.61) . is 
related to the total 57(2) conformal spin J = N + \L of the scalar operator (11.11) . Substituting 

in ( I3.2ip by its explicit expression ( I3.6p . one notices that a part of A correction to g 2 (A) 
proportional to can be absorbed into the lowest order term as follows 

q 2 (X) = [-(N + \L + |A £ (0) )(iV + \L + |Ae (0) - 1) - \L] + 0(\ 2 ) . (3.24) 

We recall that defines the one-loop correction to the anomalous dimension of scalar operators, 
Eq. ( 13.41) . Going over to higher orders in A, one finds from ( I3.2ip that Xe^ get replaced in 
( I3.24p by the multi-loop anomalous dimension ( 13.231) . Then, one deduces from ( 13.241) that the 
factor Aa Cn \x) renormalizes the "bare" conformal spin J of the scalar operator by an amount 
proportional to its anomalous dimension 

J = N+\L ^ J rcn = N + \L + | 7 (A) . (3.25) 

This result can be interpreted as follows. For a conformal operator O con f(0) = J2 n ^( n )C™(0), 
its conformal SL{2) spin J = ^(d + s) depends on its scaling dimension, d, and projection of 

lr This explains why the one-loop dilatation operator inherits the conformal symmetry of the classical La- 
grangian. Conformal anomaly affects the anomalous dimensions starting from two loops only. 
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its Lorentz spin on the light-cone, s. To order A one has d = N + L and s = N so that 
J = N + \L. To higher orders in A, the spin s is protected from perturbative corrections while 
the scaling dimension d receives the anomalous contribution 7(A) leading to J ren = N+^L+^j(X). 
One can arrive at the same conclusion from consideration of the conformal Ward identities as 
explained in details in Refs. [361 122] . In this formalism the additive correction to the conformal 
spin proportional to the anomalous dimension comes from the renormalization of the composite 
operator given by the product of the Wilson operators and the trace anomaly of the energy- 
momentum tensor in regularized gauge theory [36] . 

Let us compare (13.231) with the results of explicit diagonalization of the two-loop mixing 
matrix in the J\f = 4 SYM, Eqs. (13. ip and (12.221) . Similar to the Baxter equation, the mixing 
matrix ( I3.ip has eigenvalues with zero and nonzero quasimomentum. Although the anomalous 
dimensions of single trace operators correspond only to the former, we can perform the comparison 
for all eigenvalues. In this way, we verified that the two eigenspectra coincide for L < 5 and 
N < 20, to two loops at least. Thus, the relation (I3.23j) provides the exact solution to the spectral 
problem (13. ip for two- loop mixing matrix in M = 4 SYM. Making use of the relation (12.171) . the 
correspondence can be further extended to M = 2 theory. 

We also checked that the "odd" conserved charges q2k+i, corresponding to the paired eigen- 
values have opposite signs in agreement with the lowest order expectations, Eq. (13.71) . e.g., for 
the state with [L = 3, iV = 5] given in Table [lj one finds for the transfer matrix (13.201) 

tt =3 (u)=2u 3 ±u 2 (-±X + ^X 2 )VTTEE 

, .. ( 73 111 \ 1 7595 \2\ 1 ( 1 , 3791 \ 4894295 \2\ / —— /o 9fi \ 

+ U \~T ~ — X + "IT A ) ± U + 1584 A ~ 5018TT2 A JV1155. (3.26) 

For eigenvalues with zero quasimomentum, we summarized our results in Table [TJ At two loops, 
they agree with diagrammatic calculations of the [L = 3, iV = 2] anomalous dimension of Refs. 
[22l [20] , which is related to the BMN counterpart of the Konishi current [21] due to multiplet 
splitting, and the [L = 3, iV = 3] result of [211 [221 [26] as well as with eigenspectra of Ref. [25] 
based on algebraic construction of the dilatation operator. 



4 Discussion and conclusions 

Recently, extensive multiloop calculations in various sectors of SYM theories pointed to persis- 
tence of integrability beyond leading perturbative order [IHIITSII2QII2III221I231I211I251I2S1. In this 
paper, we continued our study of integrability properties of the two-loop dilatation operator in 
(supersymmetric) Yang- Mills theories initiated in Ref. [23]. As a case of study, we have chosen 
the sector of single-trace operators built from holomorphic scalar fields in the M = 2 and Af = 4 
SYM theories and containing an arbitrary number of covariant derivatives projected onto the 
light-cone^. To one-loop order, in both theories, the dilatation operator in this sector can be 
mapped in the multi-color limit into a Hamiltonian of the SL(2; R) Heisenberg spin chain and its 
eigenspectrum can be found by means of the Bethe Ansatz. Our goal was to understand whether 
integrability survives in high loops and if so then what are the novel features of the underlying 
spin chain. To this end, we performed an explicit two-loop calculation of the dilatation operator 
and found that the spectrum of two-loop anomalous dimensions has the same degeneracy prop- 
erties as to one loop level. We also demonstrated that, in agreement with our previous findings 

2 The twist three, spin three anomalous dimension was computed earlier in Refs. [2T1 [22l [26] . 
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[L,N] 


e (0) 




e (2) 


[3,2] 


4 


_g nnnn] 


17 na 


[3,5] 


35 
4 


18865 
1152 


1068515 
18432 


[3 81 


25 
3 

5087 
420 


455 
32 

1210695307 
49392000 


11407175 
248832 

2330723533437143 
26138246400000 


[4,2] 


5± 


-|(17±5V5) 


^(65 ±23^) 


[4,5] 


35 _i_ 

4 12 


28139 -j- 9101-v^ 
1728 + 44352 


799837 _i_ 3060649313v / 385 
13824 3688312320 




6.4113 


-8.4697 


22.4035 




9.1601 


-14.7918 


47.6639 




9.9596 


-18.2198 


62.7707 




9.8710 


-18.6154 


68.0070 


[4,8] 


12.4010 


-24.3757 


88.7702 




12.9479 


-25.7258 


93.3842 




12.9650 


-25.2831 


90.0613 




14.9761 


-30.4673 


111.0666 




16.4651 


-33.8137 


123.7385 


[5,2] 


2 
6 


3 
2 

21 

2 


37 
16 

555 
16 



Table 1: Eigenvalues of the one-, two- and three-loop dilatation operator in M = 4 SYM. 
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[23] , the two- loop anomalous dimensions in M = 2 and M = 4 theories differ from each other by 
an overall normalization factor indicating that the phenomenon is not sensitive to the conformal 
symmetry. These results lay a firm ground to the belief that the dilatation operator in the two 
theories is integrable beyond one loop. 

As a next step, we tried to uncover integrable structures behind the two-loop dilatation op- 
erator by applying the method of the Baxter Q— operator. The reason for this is the following. 
It is well known that in classical integrable models admitting the Lax representation one can 
apply the "magic recipe" [32] to perform a canonical transformation to the separated variables 
and reduce the original multi-dimensional problem to a set of one-dimensional ones. In quantum 
integrable models, the canonical transformation is replaced by a unitary transformation to the 
separated coordinates such that the multi-particle wave function is factorized into a product of 
single-particle ones. For the 5L(2;R) Heisenberg spin chain (= one-loop dilatation operator) 
the representation of the separated coordinates (SoV) has been constructed in [33]. In this rep- 
resentation, the single-particle wave function is given by the eigenvalue of the Baxter operator, 
Q(°\u), and the Schrodinger equation in the separated variables coincides with the Baxter equa- 
tion (13.101) supplemented with (13. lip and (13.91) . Going over to higher loops, we assumed that the 
spin chain describing the eigenspectrum of the multi-loop dilatation operator admits the SoV 
representation with the single-particle wave function Q(u) corrected by perturbative corrections 
(I3.16p . In the M = 4 SYM, this is in agreement with the fact that the dilatation operator for 
scalar operators with large canonical dimension L and Lorentz spin iV can be identified via the 
gauge/string correspondence with a Hamiltonian of the classical sigma-model on AdSaxS 1 back- 
ground. This sigma-model is known to be completely integrable and it admits both the Lax and 
SoV representational [38] . 

The question remains however how to construct the higher-loop Baxter Q— operator and what 
is the analog of the Baxter equation (13.101) for its eigenvalue Q{u). To answer the second part of 
this question we first verified that eigenvalues of two-loop dilatation operator calculated in the 
M = 4 theory are in agreement with the modified Bethe Ansatz equations conjectured in |34j . 
Identifying the Bethe roots as roots of Q(u) we worked out a deformed Baxter equation which 
exactly encodes the one- and two-loop spectra of anomalous dimension. Then, we demonstrated 
that the Baxter equation correctly incorporates a peculiar feature of conformal operators - the 
conformal SL(2) spin of such operators is modified in higher loops by an amount proportional to 
their anomalous dimension. From the point of view of spin chains this property implies that the 
underlying integrable model is rather unusual - the Hamiltonian of the spin chain depends on 
the total SL(2) spin which in its turn is proportional to the Hamiltonian. Still, to identify this 
spin chain one needs the explicit form of the Q— operator. To one-loop order, this operator has 
been constructed in [39|. Acting on the Wilson operators in the momentum representation O(u), 
Eq. (12.31) . the one- loop Baxter operator can be realized as an integral operator Q« (u\v) acting 
on the momentum fraction, in a close analogy with the dilatation operator V(it|i>), Eq. (12. 4p . 
The one-loop evolution kernel Y^°\u\v), (12. 6p . arises as a coefficient in the expansion of the 
kernel Qu\u\v) in the spectral parameter around u = ±|, in agreement with (13.91) . In a similar 
manner, one can translate (13.231) into the relation between the two-loop evolution kernel (12.51) 
and the two-loop Baxter operator. Simplicity of the two-loop kernel (12.101) gives us a hope that 
such operator can be constructed explicitly and the problem deserves additional studies. 

3 For an interpretation of the Q— operator in string theory see Ref. [37J. 
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Appendix 

To one- loop order, the moments of the two-particle kernel vj°\ Eq. (12.71) . are given by 

A (0)mim 2 _ q(l) jmi ™ 2 , ^1^2- @n 2 m 2 xmi +m 2 , J n l ^ n 2 \ /» I 1 



mi!m2! «2 ~ m 2 m+na |^ mi ^ m2 

with discrete step-function taking values 8 nm = 1 for n > m and vanishing otherwise. Analo- 
gously, to two loops, the moments of the two-particle irreducible kernel (12.1 5p are given by 

A (1) ™ = [S$S$ + S$ - 2(4 - WS<£] CC + "T' r m2 ni+n > 2(4 - J\f) (A.2) 

_ I C?(2) _ io(2) _ 1 fc(l) _ Q(l)\ ( cr(l) ,0C(1) _4Q« \ 1 f ™1 

Here the notation was introduced for harmonic sums 

n i n 1 2 

^ = E 7 = ^ n + X ) - ^(1) > ^ = E 71 = + X ) + y • (A-3) 

£=1 £=1 

with t/>(:e) = d In T (a?) /da; being the Euler digamma function. 

The moments of the three-particle irreducible kernel /\( 1 ) mim2m3 reac i f or n — o 



^(l)mim 2 m 3 /ii./i3- Li ni + n2 _|_ n3 "n3m3 



J/g f 1 I g(i) c(i) ^ 

^ "mini I T *J mi u ni ^m\~n\ ] 

\rii — mi / 



m0n 3 milmjlma! (ni — mi){n^ — 777.3) 



777 x <-> 7773 

while for 77 2 > they can be expressed in terms of the n 2 = moments as 

^Y(l)mirrt2m 3 _ \ ^ \ * \ ^ \ ^ I jJ n 2- j^(l)m 1 -k 1 ,m2-n 2 +j+k,m 3 -k 3 / » g\ 

nm 2 n 3 " Z^ Z^ Z^ iJioli-! (Vic — 7 — fcV «i+Jl-0,n 3 +j3 I • J 

Jl= j 3 =0 fei=0 fc 3 =0 J ^ J6 - i- 2 J >' 

where we used shorthand notations j = ji + j$ and k = k\ + k 3 . 
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